Magneto-orbital effect without spin-orbit interactions 
— noncentrosymmetric zeolite-templated carbon structure 
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A peculiar manifestation of orbital angular momentum is proposed for a zeolite-templated carbon 
system, C36H9. The structure, being a network of nanofiakes in the shape of a "pinwheel", lacks 
inversion symmetry. While the unit cell is large, the electronic structure obtained with a first- 
principles density functional theory and captured as an effective tight-binding model in terms of 
maximally-localized Wannier functions, exhibits an unusual feature that the valence band top comes 
from two chiral states having orbital magnetic momenta of ±1. The noncentrosymmetric lattice 
structure then makes the band dispersion asymmetric, as reminiscent of, but totally different from, 
spin-orbit systems. The unusual feature is predicted to imply a current-induced orbital magnetism 
when holes are doped. 
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I. INTRODUCTION 

Functional materials, as exemplified by magnetoelec- 
tronic materials, ferroelectrics, thermoelectrics, or un- 
conventional superconductors, usually rely essentially on 
strong Coulomb correlations or spin-orbit coupling in 
the system. That is why magnetic elements or heavy 
elements play a crucial role there, as in the uncon- 
ventional superconductivity involving transition-metal 
compounds^ and topological insulators involving heavy 
elements 2 -^. All the more acutely, we can pose a con- 
verse question: can we explore possibilities to realize such 
unique physical properties in light- element materials, by 
designing unconventional structures? 

Among light elements, carbon is unique and promising. 
Its covalent nature enables the atoms to have a variety 
of structures, ranging from three-dimensional diamond 
down to zero-dimensional fullerene, and the electronic 
structure dramatically changes accordingly. The scale of 
the band width goes from ~ 10 eV in graphene^ down to 
< 1 eV in solid fullerene^. Carbon nanotube becomes a 
metal or semiconductor depending on its chiralitji In- 
deed, there is a body of theoretical and experimental in- 
vestigations on functional materials made from carbon. 
For instance, Shima and one of the present authors the- 
oretically proposed that long-period (i.e., antidot-array) 
graphene structures accommodate a unique platform for 
ferromagnetism as well as a manipulation between metal- 
lic and semiconducting band structures^. 

Now, an even more intriguing avenue that may har- 
bor rich possibilities in electric properties is, in our view, 
three-dimensional periodic networks of carbon atoms. A 
recent remarkable example due to Kyotani and coworkers 
is the zeolite-templated carbon (ZTC), which is a unique 
class of carbon networks fabricated by employing a zeolite 
framework as a template (where zeolite is taken away af- 
ter fabrication)*^. A fascination with zeolite is its versa- 
tile (almost 200 known) structures^, so various ZTCs are 



expected to be possible. Another fascination is zeolites' 
open structures with large cages arranged with long pe- 
riods. In fact, regularly arranged alkali clusters in zeolite 
exhibit a variety of unusual physical properties, including 
ferromagnetis m 12 ' 13 despite the system comprising only 
nonmagnetic elements. An important theoretical picture 
that has emerged from the system is that the electronic 
structure of the alkali-metal loaded zeolites can be cap- 
tured in the "superatom" modeli^— , where the whole 
electronic band structure can be described surprisingly 
accurately as a tight-binding model of Wannier orbitals, 
which reside on each zeolite cage with a linear dimen- 
sion ~ 10 A but can have symmetries such as s and p. 
A message suggested by this is we may design versatile 
functional systems by controlling the arrangement of su- 
peratoms and/or the nature of each superatom. 

In the present work we go one step further in seek- 
ing a possibility of realizing unusual electronic structures 
from unusual geometries. We shall propose that a zeolite- 
templated carbon can indeed harbor, when its crystal 
structure is noncentrosymmetric, an unusual, asymmet- 
ric band dispersion arising from "chiral" superatoms hav- 
ing orbital magnetic momenta of ±1. So in this pa- 
per, we first obtain, with an ah initio density-functional 
calculation, the electronic structure of ZTC, for which 
the zeolite Y is used as the template. While this sys- 
tem has been experimentally studied intensively due to 
its high hydrogen storage ability*^, theoretical studies 
have been quite limited, since the crystal structure is 
extremely complicated and has not been determined ex- 
perimentally Recently, Kyotani et al. have constructed 
several possible atomic structure models^. There, one 
possibility is a kind of network of carbon nanotubes^. 
Another possibility, on which we focus here, is a network 
of nanographenes in a form of buckybowl (as opposed 
to buckyball) in a three-dimensionally regular structure. 
Here we take the simplest model, in which the build- 
ing block is C36H9 nanographene with hydrogen termi- 
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nation (called C36H9 patch hereafter) (Fig. 1(a)). Wc 
shall show that, while there are as many as 24 valence 
bands due to the complicated structure, it is possible to 
construct a simplified tight-binding model in terms of the 
maximally-localized Wannier functions for these bands, 
where the essence is a surprisingly simple three-orbital 
model on the diamond lattice with eight lattice points 
per unit cell. The three orbitals, ipi,ip2, and ^3, on each 
C36H9 patch can be seen, respectively, as p x , p y and p z 
orbitals of a superatom with the orbital angular momenta 
l z = —1, 0, 1, although the symmetry of the three orbitals 
is lower than that of the atomic p orbitals. 

Now, what is unusual is the following: the present 
crystal structure may be viewed as a network of (three- 
winged) "pinwheels" (consisting of four C36H9 patches; 
Fig. 1(b)), where each pinwheel has a threefold axis but 
lacks the inversion symmetry (hence its name). This de- 
prives the crystal structure of the inversion symmetry as 
well, which makes valence top made from chiral states, 
which arc symbolically -0g oc (p x +e ±27 ™/ 3 p J/ +e ±4,I ' l / 3 p z ), 
which belongs to the irreducible representation of E (sec 
below). The two chiral bands are shown to have differ- 
ent band dispersions, which reminds us of spin-orbit split 
bands in GaAs but the cause of the splitting in the for- 
mer is totally different. We finally predict that we can 
induce an electric-ficld-induccd orbital magnetism if wc 
dope carriers into the system and apply an external elec- 
tric field. Although the asymmetric band dispersion has 
nothing to do with spin-orbit interaction, the unusual 
feature implies, in its own right, that a current-induced 
orbital magnetism is expected by hole doping. 



II. METHOD 

To obtain the electronic properties, we have performed 
a first-principles calculation 1 ^ within the local-density 
approximation (LDA) 20 i 21 . Ultrasoft pscudopotcntials 
are used to describe the electron-ion interaction^ 2 -. The 
valence wave functions and charge densities are expanded 
by a plane-wave basis with cutoff energies of 30 Ry and 
150 Ry, respectively. Integration over Brillouin zone is 
carried out with 2 x 2 x 2 k points. We have employed the 
wannier90 code^ 3 . to compute the maximally-localized 
Wannier orbitals. 



III. RESULTS 

Figure [T] shows the geometry of C^flg. The struc- 
ture is theoretically optimized with the experimental lat- 
tice constant of a = 24.07 A. Since the zeolite Y tem- 
plate has a diamond-type structure, the resulting C36Hg 
is diamond-like as well. However, unit cell of the whole 
structure is simple-cubic rather than face-centered cubic, 
and contains eight hydrocarbon C36H9 patches, whose 
normal directions are all different, i.e., (±1,±1,±1). Each 
C36H9 patch, connected to three nearest-neighbor ones as 




FIG. 1. (Color online) (a) Crystal structure, theoreti- 
cally optimized here, of the zeolite-templated carbon. Solid 
line represents a unit cell, which contains eight hydrocarbon 
C36H9 patches (with 288 carbon (darker spheres) and 72 hy- 
drogen atoms (lighter)) with different orientations, (b) A 
C36H9 patch with its three nearest-neighbor ones, which form 
a three-winged "pinwheel" with a threefold axis. 
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FIG. 2. LDA band structure of C36H9. Origin of energy is 
set to be the top of the valence band. Inset shows the first 
Brillouin zone. 

shown in Fig.JIJb), form a three- winged "pinwheel" with 
a threefold axis. This is how the crystal structure has a 
high symmetry (space group: P4i32), but lacks inversion 
symmetry. 

The obtained electronic band structure of C36H9 in 
Fig. [2] shows that C^flg is an insulator with a band 
gap of 0.84 eV. Since the three-dimensional network com- 
prises nanoflakes, the band structure resembles those of 
molecular crystals where weakly coupled molecular or- 
bitals result in a narrow band widths. The valence band 
top (— 1 < E < in cV) comprises 24 bands that are 
isolated from lower ones with an energy gap. Interest- 
ingly, the top and bottom of this region are delineated 
by almost flat bands. The band widths of these nearly 
flat bands are around 0.05 eV, which is much smaller 
even when compared with that of ti M -band in fullerene 
compounds (~ 0.3 - 0.5 eV). Naively, one might expect 
magnetic or other instabilities with hole doping for these 
flat bands, but we have not found any magnetic orders 
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such as fcrromagnctism or antiferromagnetism within the 
local spin-density approximation (LSDA). 

Let us then characterize these bands in terms of 
maximally-localized Wannicr functions. Figures [3][a)- 
(c) visualize the maximally-localized Wannier orbitals for 
the 24 bands. On each "pinwheel" reside three, equiva- 
lent Wannier orbitals that spread, respectively, along the 
three "wings" (i.e., along the three next-nearest direc- 
tions). These three orbitals, which we call (ipi,ip2,ip3)y 
conform to the threefold symmetry of the pinwheel, 
and, in the superatom language, roughly correspond to 
(p x ,Py,Pz) orbitals. Thus we may regard the system as 
a "p-electron superatom system" on the diamond lattice. 
However, crucial differences between the present super- 
atom system and usual diamond are: (i) there are no 
s-orbitals, and (ii) p-orbitals do not have full symmetries 
that the atomic p-orbitals have and are symmetric only 
under a three-fold axis rotation. This can be seen from 
the fact that each p-orbital has its center of wave func- 
tion shifted from the center of the superatom (i.e., the 
center of the pinwheel) unlike the atomic p-orbitals. 

Due to the symmetry of the present p-orbitals, these 
valence bands should be characterized by an A state and 
two E states given as 

= ^ W>i + V>2 + V>s) , (1) 
^ = ± + e±^ 2 + e ±4 "/^ 3 ) , (2) 

where ipx Qz = state) is the symmetric combination of 
the three Wannier functions, while ip^ (l z = ±1 states) 
have finite angular momenta associated with rotation 
around the threefold axis. With these states it is found 
that, the doubly-degenerate valence-top wavefunctions, 
for example out of the 24 bands, comprise entirely from 
ip^ at the r point from the symmetry analysis 2 ^. The 
calculated weights of ipA and V'e are displayed in Figs. 
IDJa) and (b), respectively. Interestingly, the nearly flat 
bands along the top and bottom of the 24 bands origi- 
nate from the finite-angular-momentum tp^: states, while 
the dispersive bands from tpx- 

Figure [5] illustrates the tight-binding model for the 
24 bands obtained in terms of the maximally-localized 
Wannier orbits (wings of a pinwheel in the figure). 
Since the maximally-localized Wannier functions are ba- 
sically localized on each patch, it is natural that the 
3 x 8 = 24 bands are all accurately described by them. 
There are three transfer integrals, t\ = —0.255 eV, 
t 2 = 0.116 eV, and £3 = 0.022 cV, between nearest- 
neighbor patches, while the largest transfer integral be- 
tween second-neighbor patches is £4 = —0.092 eV. These 
exhaust the transfer integrals whose magnitudes exceed 
0.01 eV. Although the unit cell of the ZTC is huge and 
complicated (containing more than 300 atoms), these 
four transfer integrals suffice to reproduce the dispersions 
of the 24 valence bands around the Fermi level. 



IV. SEPARATED DISPERSION AND 
CURRENT-INDUCED ORBITAL MAGNETISM 

Now we come to the key result in the present work 
in Fig[6j which depicts the structure of the topmost five 
valence bands, with the main character of each band indi- 
cated in terms of tpA , V'e ' an d V'e ■ Remarkably, the two 
topmost bands are asymmetric with respect to T point. 
They arise from the absence of inversion symmetry in the 
crystal structure, and transform with each other by the 
time-reversal. 

The feature may seem reminiscent of the noncen- 
trosymmctric systems with strong spin-orbit coupling, 
but the cause of the splitting is totally different here. 
A conventional spin-orbit model is Rashba's 2 ^: in a two- 
dimensional electron gas with a strong spin-orbit cou- 
pling, the spin is no longer a good quantum number, 
but the bands are still degenerated due to the Kramers 
degeneracy. If we apply an electric field perpendicular 
to the plane, the degeneracy is lifted, and a vortex-like 
spin texture is generated. Even in the absence of electric 
fields, various spin textures can emerge if the system is 
inversion asymmetri o 26 i 27 . In fact, the original proposal 
by Rashba was made for noncentrosymmctric wurtzite 
semiconductors 2 ^—. Such systems harbor a possibility 
of manipulation of the spin degree of freedom by applying 
an external electric fiel d 31 i 32 . In materials belonging to 
the gyrotropic point group, charge current is generically 
accompanied by a nonzero spin polarization 3 - 3 -. 

By contrast, in the present ZTC we can expect a 
current-induced orbital magnetism: If we dope holes to 
extract the peculiar property of the valence top and ap- 
ply an external electric field, the charge carriers drift in 
the direction of the applied field, so that the popula- 
tion of the l z = 1 states will exceed that of l z = — 1 as 
in the spin-orbit case 3 ^, since the band structure of the 
two orbital-magnetic states is asymmetric around the F 
point as shown in Fig. |Ef 5 . For example, if we dope 
0.05 holes per unit cell and apply an external electric 
field in (1,1,1) direction so that the electrons acquire the 
quasimomentum Ak = 2 x 10~ 2 in the unit of 2tt /a, the 
resulting orbital magnetic moments, in units of 10 -4 /ib 
per unit cell, are = 0.47, (fir_x 1 i\) = 

-0.50, = (M(i.-i,i)> = = 0-40, 

and _ Xjl )) = (M(-i,x -1)) = = -0.39, 

where ((J-( a ,b,c)) denotes the orbital magnetic moment 
contributed from C36H9 patch whose normal direction 
is (a,b,c). Since the quantization axis differs from one 
patch to another, these moments do not cancel with each 
other, and the total orbital magnetic moment in (1,1,1) 
direction is (fx) = T, a ,b,c=±i(l J >(a,b,c))(a + b + c)/3 = 
1.8 x 10~ 4 /iB- When the electric field and/or the num- 
ber of holes are varied, the magnetic moment increases 
roughly proportional with the field and with the number 
of holes as shown in Fig. |6tc) . 

The current-induced orbital magnetism in the present 
mechanism exploits the fact that electrons are confined 
on the curved surface of the C36H9 patch. Namely, in an 
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FIG. 3. (Color online) Maximally-localized Wannier functions, 4>i, -02, and ip3, 
band. The surfaces with different hues represent positive (negative) isosurfaces. 



for the 24 bands at the top of the valence 





FIG. 4. Structure of 24 topmost valence bands with the line 
width representing the weights of (a) ipA state, or (b) ip^ . 



FIG. 5. (Color online) The C36H9 represented by the ob- 
tained tight-binding model. Each black sphere represents a 
C36H9 patch, while wings of a pinwheel illustrate Wannier or- 
bitals. The transfer integrals are: ti = —0.255 eV, £2 = 0.116 
eV, t 3 = 0.022 eV, and U = -0.092 eV. 



the usual exploitation of magnetoelectric effects in solids 
for spintronic o 36 ' 37 , where appreciable spin-orbit interac- 
tion is required to couple the motion of charge with the 
spin degree of freedom. The present mechanism for the 
current-induced orbital magnetism requires no spin-orbit 
coupling, so that we may expect that this mechanism will 
open a new avenue for light-clement functional materi- 
als. Further systematic studies on electronic structures 
of ZTC systems is an interesting future problem. 



electric field along, say, (1,1,1), electrons on a pinwheel 
_L (1,1,1) cannot move along that direction, but edge 
currents are generated instead. Due to inversion asym- 
metry, the right-handed and left-handed edge currents do 
not cancel with each other, which will induce an orbital 
magnetic moment. Thanks to the small width of the rel- 
evant bands, a relatively small electric field should suffice 
to induce such an orbital magnetism. This contrasts with 



V. CONCLUSION 

We have found that the zeolite-templated carbon sys- 
tem, C36Hg, has a unique electronic structure around 
the top of the valence band, where orbital character is 
asymmetric around Y point due to a combined effect of 
the chiral nature in the orbital character and the lack 
of inversion symmetry. As a consequence of this char- 
acteristic band structure, a current-induced orbital mag- 
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netism is prcdected when holes are doped. An orbital 
magnetic texture, which resembles spin textures in rel- 
ativistic noncentrosymmetric systems, is then expected 
to appear despite the absence of spin-orbit coupling, and 
can be understood in terms of the simple tight-binding 
model. 
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FIG. 6. (Color online) (a) A blow-up of Fig(4] where bands 
originating from ipA (black dotted line), ip£ (red solid), and 
(blue dashed) on a patch normal to (1,1,1) are displayed, 
respectively, along T - R line, (b) Two topmost valence bands 
plotted on a two-dimensional k— space. The two bands anti- 
cross with each other along the (1,1,1) direction, and arrows 
represent the orbital moment, l z . (c) The total orbital mag- 
netic moment per unit cell as a function of the electric field 
with the number of holes, 8 = 0.01, 0.02, 0.03, 0.04 and 0.05 
(from bottom to top). The applied electric field is here rep- 
resented by Ak, which is the crystal momentum acquired by 
an electron due to the field in (1,1,1) direction. 



